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Abstract 

Grand unified theories based on large groups (with rank > 6) are a natural context for dark 
matter models. They contain Standard-Model-singlet fermions that could be dark matter candi¬ 
dates, and can contain new non-abelian interactions whose sphalerons convert baryons, leptons, and 
dark matter into each other, “cogenerating” a dark matter asymmetry comparable to the baryon 
asymmetry. In this paper it is shown that the same non-abelian interactions can “pre-annihilate” 
the symmetric component of heavy dark matter particles %, which then decay late into light stable 
dark matter particles ( that inherit their asymmetry. We derive cosmological constraints on the 
parameters of such models. The mass of x must be < 3000 TeV and their decays must happen 
when 2 x 10“^ < Tdec/m^ < 10“^. It is shown that such decays can come from d = 5 operators 
with coefficients of order I/Mqut or 1/Mp£. We present a simple realization of our model based 
on the group SU{7). 


1 Introduction 

In recent papers it has been argued that grand unification based on groups larger than SU{5) has 
several features that make it a natural context for theories of dark matter [H [2]- (1) The fermion 
multiplets of such unified groups typically contain Standard Model singlets, which could play the role 
of dark matter. (2) Unified theories can have accidental global symmetries (a well-known example 
is n — L in minimal SU{5)), and one of these could be the quantum number X that stabilizes dark 
matter. And (3), unification groups of rank 6 or larger can contain additional non-abelian interactions, 
which can play a role in generating dark matter. One way this can happen is through the sphaleron 
processes [3] of the new non-abelian interaction, which can convert baryons and leptons into dark 
matter particles, thus giving a scenario in which dark matter and baryonic matter are “co-generated.” 


Many papers have explored the idea that sphalerons are involved in the genesis of dark matter. 
Usually, these scenarios involve only the sphalerons of the electroweak interactions 015]. However, 
as emphasized in [6] (and later in [2]), the idea that the sphaleron processes of a new non-Standard 
Model interaction do the co-generation has several attractive features. First, if electroweak sphalerons 
produced the dark matter asymmetry, the dark matter particles have to be chiral under electroweak 
SU{2), which makes them more visible both directly and through their loop effects, severely constrain¬ 
ing such models. By contrast, to be produced by the sphalerons of a new interaction, the dark matter 
particles need only be chiral under the new interaction and can be vectorlike under the electroweak 
SU{2), which allows them to be heavy and to have negligible loop effects on Standard Model predic¬ 
tions. Such models are therefore less severely constrained and much easier to build. Second, having 
two kinds of sphalerons — those of the electroweak SU{2) and those of a new gauge interaction — 
means that there are two sphaleron equilibrium conditions, which is enough to determine the ratios 
of B, L, and X. That means (under the assumption that the dark matter is almost purely asym¬ 
metric, i.e. <C tidm) that the number density of dark matter particles today can be calculated, 

giving a definite (though of course model-dependent) prediction for the mass of the dark matter par¬ 
ticle. As sphaleron processes typically make B, L, and X comparable, and because and ^dm are 
comparable, one typically finds that the mass of the dark mass particle is GeV-scale . 

To have the dark matter be asymmetric [HliaElEllTllgKio], there must have been some processes 
in the early universe that annihilated the “symmetric component” of the dark matter. That is, most 
of the anti-dark matter particles must have annihilated leaving just the asymmetry (as was the case 
with baryons). This is a non-trivial issue if the dark matter particles have masses of order a GeV, 
as then annihilation via electroweak interactions is not efficient enough and would have left a relic 
abundance of anti-dark matter particles that was far too large. One way to solve this problem [2] is to 
posit the existence of a light field (typically it must have mass less than the weak scale) that mediates 
an interaction by which dark matter particles and their anti-particles can annihilate, as well as some 
new massless (or extremely light) particles into which they can annihilate. Thus the existence of an 
entire sector of light particles, besides the dark matter particles themselves, would be postulated ad 
hoc. 

In this paper, we point out another way: The same new non-abelian gauge interactions whose 
sphalerons co-generate dark matter in the first place, are also able to annihilate the symmetric com¬ 
ponent of the dark matter. This leads to a very economical, tightly integrated scenario, in which 
everything comes from unification in larger groups: the dark matter fermions themselves, and the new 
interaction that generates them and annihilates the symmetric component. 

At first, it might seem that the new non-abelian gauge interactions whose sphaleron processes 
generate the dark matter asymmetry cannot do the annihilation efficiently enough, any more than 
the weak interactions could, for the scale characterizing their breaking must be even larger than the 
weak scale. However, as we show in this paper, they can do it through a mechanism that we call “pre¬ 
annihilation.” This mechanism was proposed in the context of other scenarios by Kitano and Low m 
and by Kang and Li [13j . (It can also be seen as a limiting case of the “exodus” mechanism of Unwin 
|14j.l Because the scenario we are describing is tightly integrated and to a large extent determined by 
the grand unified structure, one can use cosmological observables to constrain it. 

The idea of pre-annihilation is that an asymmetry in X is first produced among some heavy 
metastable dark particles, which are able to annihilate efficiently enough to make negligible the number 
of their anti-particles. These heavy dark particles then decay late, out of equilibrium, to stable dark 
particles that “inherit” their asymmetry and form the dark matter that is observed today. (We shall 
call the heavy dark particles x the light stable dark particles C-) We shall show that the lifetime 


of these late decays must fall within a certain window. A nice feature of the scenario we are discussing 
is that these decays can result from effective dim-5 operators that would give the right lifetime if they 
were suppressed by 0{1/Mpi) or 0{1/Mgut)- 

The basic cosmological scenario we shall discuss has the following stages. When the temperature 
is above the mass scale of the heavy dark particles, (which is much larger than the weak scale but 
less than 3 x 10® GeV), the sphaleron processes of the new non-Standard Model gauge interaction, 
whose gauge group we shall call G*, produce an X asymmetry comparable to the baryon asymmetry. 
This asymmetry is in the form of an excess of y over their anti-particles x^- As the temperature drops 
below m^, the x x^ annihilate by means of the G* gauge interaction, leaving only the excess x 
due to the X asymmetry. At a later epoch (when T is between and 2 x each heavy 

dark particle y decays via interactions suppressed by 1/Mgut or \/Mpi into one light dark matter 
particle ^ plus some Standard Model quarks and leptons that equilibrate. The C are singlets under all 
gauge interactions — except possibly ones broken at Mqut — and are never in thermal equilibrium. 
(At some very early time, their density was driven to exponentially small values by inflation, and the 
only such particles present today are those caused by the decay of the heavy dark particles.) Thus 
the stable light dark matter particles inherit the asymmetry produced at very large scales. Because 
they have never been in equilibrium the dark matter particles are produced with a mean momentum 
significantly higher than the photon temperature. Depending on the parameters of the model, this 
scenario can give either cold or warm dark matter. 

As noted above, some of the ingredients of the scenario we are proposing and studying here have 
been considered elsewhere, such as the use of non-Standard Model sphalerons to co-generate dark 
matter PIS], and the pre-annihilation of the symmetric component of dark matter [EllISKIl]. What 
is new here is that we combine these elements into a relatively simple scenario and show that it emerges 
naturally in the context of grand unification based on groups of rank 6 or greater. And within models 
of the type we are proposing, we derive significant new limits coming from several considerations. 
These include (a) the condition that the decays of the heavy unstable dark matter not be so slow as to 
make the resulting stable dark matter hot, and (b) that the decays of the heavy unstable dark matter 
not be so rapid as to regenerate a large symmetric component of dark matter. These limits give a 
relatively narrow window of decay lifetimes of the heavy dark matter, which we show turns out to be 
nicely consistent with decay via GUT- or Planck-suppressed operators that one might expect to exist 
in realistic models. 

All of this we illustrate in a specific grand unified model based on SU{7). In that model, we show 
how the various critical ingredients tend naturally to exist in grand unified models with large gauge 
groups. The model illustrates, in particular, how the dark matter and ordinary matter particles are 
unified in GUT multiplets; how the G* interactions are unified with the Standard Model interactions; 
how in a typical case there are quarks, leptons and dark matter particles that transform under G*, as 
needed for the co-generation mechanism; and how the quantum number X can arise as an accidental 
symmetry in a grand unified model. 

We shall first discuss the cosmology of the scenario in more detail in section 2, and then the details 
of the illustrative SU{7) model in section 3. 


2 The details of pre-annihilation 


The new gauge interaction based on the group G* is broken at a scale M* by the VEV of a Higgs 
field that we shall call H. It is the sphalerons of the G* interaction that will violate X and relate the 
X asymmetry to the B and L asymmetries. It is also the G* gauge interaction that will mediate the 
annihilation of the heavy dark particles (y) with their anti-particles The heavy dark particles 

(x) which have X ^ 0, are neutral under the Standard Model gauge interactions, but transform in 
a chiral way under G*. (It is because of this chirality, of course, that the G* sphalerons violate X.) 
The particles x and obtain masses of order M* from the VEV of (H). We can therefore denote the 
scale of the breaking of G* by (instead of M*) and shall do so from now on. Some of the Standard 
Model quarks and leptons will also transform under G*, so that G* sphalerons also violate B and L. 
(A simple model having all these features will be presented in section 3.) 

The light dark particles (C) and anti-particles (C'^) are singlets under SU (3)c x SU (2) l x U{1)y x G*. 
They couple to the other particles of the model through dimension-5 operators suppressed by 1/Mcut 
or 1/Mp£. 

We assume that at very early times inflation wiped out B, L, and X. After (or during) reheating 
an asymmetry of one or more of these quantum numbers was generated. It will not matter for our 
discussion which asymmetry was generated or what process generated it. We assume that no significant 
number of C and particles appeared during reheating, because they have no coupling to the inflaton 
and only couple to other light particles through GUT-scale-suppressed operators. (This will put an 
upper limit on the reheat temperature that depends on the GUT scale, as will be discussed later.) 

When T was large compared to m^, the G* sphaleron processes were happening and produced an 
equilibrium among B, X, and L, making these asymmetries comparable. (Of course, as is common, 
we use the term “sphaleron process” loosely to mean any anomalous G* interaction, not just the ones 
involving tunneling via the sphaleron conhguration.) Let us define ep = (np —n-g)/s, ep = (n^ —n-jj)/s, 
and ex = {nx — n-j^)/s. In order to have ^dm ~ 5Hp, one requires that ex ~ (5mp/m(^)ep, where 
nip is the proton mass. 

When T dropped below m^, the number of x and x'^ became Boltzmann suppressed as they 
annihilated into lighter states. These annihilations would have happened primarily through G* gauge 
interactions (and partly through Q exchange). (Therefore, they would not have annihilated into the 
light dark particles C + C'^-) The relic abundance of an asymmetric particle can be expressed in terms 
of the asymmetry parameter r, which is dehned to be 


n: 


r = 




( 1 ) 


Eor an annihilation cross section scaling as ao{m^/T)~ 
Plo] (see also earlier calculations in Refs. HSj): 


the value of r at late times, r^o, is given by 


^ OO 



^x^V9*f \ 
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( 2 ) 


where A = 0.264 Mpum-^ao (with Mpi denoting the Planck mass), g^f is the effective number of 
relativistic degrees of freedom when the annihilations freeze out at temperature Tj, and Xf = m^/Tf. 
(We will use throughout to denote relativistic degrees of freedom, with the subscript labeled to 
denote the temperature at which this quantity is evaluated in each case). Taking the annihilation in 



( 3 ) 


our model to be s-wave (n = 0) with ao = a^^lrr? ^ Eq. ([2]) simplifies to 


exp 


0.264 e^a,Mpe^ \ 

m^Xf I 


Our goal is to achieve ~ 10“^^ (the same order of magnitude as for baryons). Somewhat 
arbitrarily, we will take the dividing line between symmetric (roo = 1 ) and asymmetric (roo ^ 1 ) 
freeze-out to be the point at which the argument in the exponential of Eq. m is — 1 ; asymmetric 
freeze-out will occur when the absolute value of this argument is > 1. For the model parameters 
considered here, y/g*f ~ 20 and xj ~ 20. (As noted in Ref. [9], xj- is little changed from its symmetric 
value unless Too ^ !)■ Assuming that a* ~ 10“^, we hnd that the y freeze-out will be asymmetric as 
long as < 3 x 10® GeV = 3000 TeV. When this limit is satisfied, to a good approximation the only 
heavy dark particles x that remain after freeze-out are those that existed due to the X asymmetry. 

As mentioned, we assume each x decays via a dimension-5 GUT suppressed operator into a light 
dark particle C plus Standard Model particles. Let us call the coefficient of this dimension-5 operator 
c^/Mgut, where Mqut is the scale where the Standard Model gauge couplings unify, and C 5 is a 
dimensionless number. The x decays will happen when the temperature of the universe is Trfeo given 
by 
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(4) 


So far, we have been giving a fairly standard analysis of asymmetric freeze out. We now show that 
there are two considerations that yield both an upper bound and a lower bound on Tdec that must 
be satisfied to achieve a satisfactory asymmetric dark matter scenario. The upper bound comes from 
considering the x'^ that decay into when T ~ m^, when the x'^ are not Boltzmann suppressed. Of 
course, only a small fraction do, since these decays have a lifetime much larger than the age of the 
universe when T ~ m^. But this fraction must be smaller than ex = ~ 10“^® if one is not 

to produce a number of that is comparable to the C asymmetry. This gives a lower bound on the 
lifetime of the decay, and thus an upper bound on Tdec- One has 


(5) 

\ S J tdec \ S J tdec g*m^ 

where is the age of the universe when T ~ m^^, and tdec is the age of the universe when T = Tdec- 
Demanding that this be less than ex — gives 


tm^ 1 iy/ 1 Buip _ Tdec (5mp\ g*m^ 1/2 i n—4 

— 7 . r z \ \ I I * \^/ 

tdec 9*m^ \y/9*de(X^^(,/Mpi) m^ y Ul^ j gJdec 

We have neglected two additional effects that will tend to suppress the decays of the . At T ~ 
the decay rate will be Lorentz suppressed, and inverse decays will be energetically allowed. Both of 
these will decrease the number of decay-produced so our limit is a conservative one. 

The lower bound on Tdec comes from requiring that the dark matter streaming length be small 
enough to allow for the formation of the observed large-scale structure. When the dark matter particles 


















C are produced by decays of the y, they have a momentum that is of order m,^/3 (assuming three-body 
decays). This momentum redshifts as the inverse of the scale factor, so a lower T^ec allows for less 
redshift prior to the epoch of structure formation, resulting in a larger velocity and corresponding 
larger free-streaming length. 


The comoving free-streaming length A ps up to the time tgq when the matter and radiation densities 
become equal and density perturbations can begin to grow, is given by m 


^FS 



(7) 


where R is the scale factor, and the zero subscript refers to present-day quantities. Using the ap¬ 
proximation that the velocity is u = 1 when is relativistic and scales as 1/R when it becomes 
nonrelativistic at tjsfR, we obtain m 
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The momentum of C at scale factor R is given hy p = {m^/2>){Rdec/R)-, and Q becomes nonrela¬ 
tivistic when p ~ m^, so we have 

1/3 


TnR — Tdec 


3m^ \ / g^dec \ 
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(9) 


where we have used the fact that g^:{RT)^ is constant [16]. In the radiation-dominated era, t = 
O.SOg'* so we obtain 

/ \ 2 


tNR = O.OSMpf 
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Further, Rq/R^ji in Eq. ([7|) can be written in terms of Tdec- 

Ro _ 3mc fTdec\ (9*decy^^ /liy/3 
Rnr m^KToJUs/Aj U; ’ 


( 10 ) 


( 11 ) 


where the two factors on the right-hand side include the heating from the annihilation of all particles 
heavier than electrons, and the heating due to electron-positron annihilation, respectively. 

Assembling all of these factors and using the measured values for teq and Tq, the free-streaming 
length can be written as 
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A reasonable order-of-magnitude upper bound on Xps to allow structure formation Xps ^ 100 
kpc. Using this limit and appropriate values for the degrees of freedom, we obtain: 

^ > 2 X 10-" ( . (13) 

Roughly speaking, values of Tdec/^n^ no more than an order of magnitude above this bound will 
produce C that behaves as warm dark matter, while larger values of Tdec/Fix yield cold dark 

matter. 





















Clearly, there is a broad range between the limits given in Eqs. ( 6 ) and (|13ll . Substituting Eq. (4) 
into those inequalities gives a range for Mgut / C5 


2 X 10^® GeV 



m-v 


3 X 106 GeV 


^ Mgut 
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< 8 X 10^® GeV 
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\3 X 106 GeVy ■ 
(14) 


We have assumed that a negligible number of are generated by scattering of thermal particles 
after reheating. In particular, we require that n^c/s < ex = {5'mp/mQ)eB- This gives the following 
upper limit on the reheat temperature: 


\ S J Trh ^GUT 


^ ^ fr. / ^ 

— < [5mp/m^)eB 

yy^rh-^Y-ji y^rh 


{Mgut/c^Y ( 5mp\ 

3 X 1025GeV V y ■ 


(15) 


All of the bounds we have derived in this section are based on the assumption that there is a strong 
asymmetry in the dark matter, i.e. <C udm- We regard this as an interesting limit. Moreover, in 

models where the dark matter asymmetry is generated by sphalerons of a new non-abelian interaction, 
assuming that limit typically allows a prediction for the mass of the dark matter particle, as noted in 
[ 6 ]. But it is interesting also to consider the more general case, where there is a weaker asymmetry 
in dark matter. This would allow smaller masses for the light stable dark matter particle C- It would 
also allow M* (and m^) to be larger, as less efficient pre-annihilation is allowed. It would also relax 
the upper bound on Tfiec/'mx- We leave a detailed analysis of the more general case to future work 

m- 


3 An SU{7) example 

We now present a simple SU{7) grand unified model that illustrates the ideas discussed above. As 
noted in the Introduction, grand unification based on large groups has several features that make it a 
natural context for models of dark matter. (1) The large unified fermion multiplets typically contain 
SM singlets that could be dark matter candidates, so that such models would unify dark matter with 
quarks and leptons. (2) The symmetry that stabilizes the dark matter can arise as an accidental 
global symmetry in a manner analogous to the way U{1)b-l arises as an accidental symmetry in 
minimal SU{5). And (3), the large unified groups (if they have rank > 6 ) contain extra non-abelian 
subgroups that can play the role of G* in co-generating and pre-annihilating dark matter. That is, 
the G* interactions can be unified with the Standard Model interactions. All of these features are 
realized in the SU{7) model we present in this section. 

The group SU{7) contains the subgroups SU{7) D SU{5) x S'[/(2)* x U{1)tj V Gsm x <S'C/(2)* x 
U{l)Tr, where Gsm is the Standard Model gauge group, and where the generator T 7 of the extra abelian 
factor can be written (in the fundamental representation of SU{7)) as Tj = diag{j^ |, |, ——f)- 

We assume that SU{7) breaks to Gsm x SU{2)^ x U{1)tj at superlarge scales. 

The fermions of one family consist of 21 -|- 3 x 7 , which we denote ~ 

The indices etc. are SU{7) indices, while k is just a label distinguishing two of the anti¬ 

fundamental fermion multiplets. This is the smallest possible anomaly-free set that gives one SM 













family. In addition, each family is assumed to have three 5?7(7)-singlet fermions that will be denoted 
simply where m = 1, 2,3 is a label distinguishing them. Besides the Higgs fields that break 

SU{7) at superlarge scales, there are Higgs multiplets /i'^, and A: = 1, 2. These will break 

the electroweak SU{2)l and the 517(2)* and give mass to the fermions, as will be seen. (Here too the 
index k is just a label distinguishing the two 17^4 multiplets.) 

We shall use lower-case Greek letters for indices of the 517(5) subgroup of 517(7), so a,/3,7,... = 
1,...,5. We shall use capital Latin letters from the middle of the alphabet for indices of the 517(2)* 
subgroup of 517(7), so /, J,... = 6, 7. We suppress all family indices. The fermion multiplets decompose 
under 517(5) x 517(2)* as shown in Table I. In that table, the 517(7) fermion multiplets of one family 
are listed across the top row. In the column under each 517(7) multiplet are given the 517(5) x 517(2)* 
multiplets contained in it. The superscripts 0, ±1, ±2 on the 517(5) x 517(2)* multiplets are the values 
of the global charge X, which will be discussed shortly. The three rows denote three types of fermion; 
“SM” stands for Standard Model fermions, which have X = 0; “VL” stands for vectorlike fermions, 
which have X = ±1; and “DM” stands for dark matter fermions, which have X = ±2. 

Table I: The decomposition under 517(5) x 517(2)* of the fermion multiplets of a family. The 
superscripts 0, ±1, ±2 are the values of X. 



^[AB] ^ 21 

i^A = 7 

=2x7 

^/,M =3x1 

SM 

= (10,1)0 

= (5,1)0 



VL 

cT 
^ IS 
II 

■07 

= (l,2)+i 

/ [k) 

= 2 X (5, i)+i 


DM 

^[iJ] 

= (1,1)-" 


/ (^) 

Wi 

= 2 X (l,2)+2 

= 3x (l,l)-2 


Note that, as is required for co-generation of dark matter and ordinary matter by 517(2)* sphalerons, 
there are dark matter particles that transform under 517(2)* (namely those in (1,2)"*“^) as well as 
quarks and leptons that do so (namely those in (5,2)“^). This is not special to 517(7), but would be 
a generic feature of models in which 517(2)* is unified with Gsm in a simple group. 

The Higgs multiplets are assumed to have the vacuum expectation values (VEVs) (/i^®'^) = ~ 

M\y, {h' 2 ) = Vd ^ X[\y, and (Dg^^) ~ ~ Tf*. (The indices of the SU{2)l electroweak subgroup 

of 517(5) are a = 1,2, while the color indices are a = 3,4,5.) The VEV Vu will give mass to the up 
quarks; the VEV Vd will give mass to the down quarks and charged leptons; and the VEVs of D) ^ 
will give mass to the vectorlike fermions and to the dark matter, as well as breaking SU (2)*. To align 
these VEVs there must be a coupling connecting these Higgs multiplets. A coupling sufficient to do 
this is . 

The Yukawa couplings are as follows 




Assumed not to be present are couplings similar to those in Eq. (16), except with -H- It is 

easy to check that there is an accidental 17(1) x 17(1) global symmetry consistent with all the Yukawa 
couplings in Eq. (16) and the Higgs self-coupling We shall be interested in a 




subgroup of this, which we will call U{1)r x K, where K is a Z 2 symmetry. The R charges and K 
parities of the various fermion and Higgs multiplets defined so far are given by 






(17) 


When the = (1,2) get non-zero VEVs of order M*, the group S'f7(2)* is spontaneously broken, 
and the group U{l)Tf x U{ 1 )r is broken down to a global symmetry that we shall call U{l)x, whose 
generator is 


X = T 7 + ^R. (18) 

One easily sees from Eqs. (17) and (18) and the dehnition T 7 = diag{j, |, |, |, |, ——|) that the 
SU{5) X SU{2)^: multiplets in Table I have the X values shown there. One also sees that the Higgs 
VEVs (h^®'^), (/ 12 ) and (Hg^^) ~ {^ 7 ^'^) do not break ?7(l)x- 

The masses of the Standard Model up quarks come from the coupling (' 0 ["^?/)'>''^)(/i^®d^, Those of the 
Standard Model down quarks and charged leptons come from (^[“^l^/iQ,)(/i 2 ), with a being a color index 
(3,4,5) for the down quarks and a weak index ( 1 ) for the charged leptons. (Of course, in a realistic 
model of quark and lepton masses there must always be other terms to avoid the “bad” unification 
relation that the charged leptons and down quarks have equal mass at the GUT scale. But there are 
several ways of adding such terms, and this does not affect anything else we say here.) The heavy 
vectorlike fermions = (5,2)“^ and = 2 x (5, l)"''^ (see the row VL in Table I)) obtain 0(M*) 
masses from ^). The Yukawa couplings so far mentioned leave the vector like fermions 

V’/ = (1,2)'’'^ massless. (The coupling ('01^^?/)/) (/ 12 ) mixes them with the heavy neutrinos, but leaves 
the resulting mixture massless.) The masslessness of these fermions would not be a phenomenological 
problem, but they can be given Dirac masses by introducing some gauge singlet fermions into the 
model. 

The heavy dark matter particles get mass of 0(M*) from and 

The (of which there are four per family, as A: = 1,2 and I = 6 , 7) are what we called x in earlier 

sections. The (of which there is one per family) and the (of which there are three per 

family) are what we called x^- 

In order for the pre-annihilation mechanism discussed in previous sections to work, the x and x^ 
must be able to decay to light stable dark matter particles through dim-5 operators that are suppressed 
by 0{1IMqut) or 0(1/Mp£). One such operator is C)hy 2 ijTh'*‘^/M gut■ Since is one of the 
X^, as noted in the previous paragraph, and and h '2 are just the weak Higgs doublets and 

hd whose VEVs do the electroweak breaking, we may write this operator as {x'^C)h*h'*/M gut■ The 
5t/(7)-singlet fermion field C can be assigned R^ = 14+ (implying by Eq. (18) that it has V = 2 as 
the other dark matter particles V’/ ^ do), in which case this dim-5 operator conserves R, K and X. 
In order to give a Dirac mass to C) let us introduce an SU(7)-singlet fermion that has R^ = —14+ 
and therefore X = —2. (Notice that K parity prevents the C from having a mass term with and 

from mixing with 

The Dirac fermion made up of ( and must have mass of order 1 GeV. The smallness of this mass 
can be explained in a way that is technically natural m by assuming it arises from the dynamical 
breaking of a symmetry. Let us suppose, therefore that there a Z 2 symmetry, which we will call P, 


under which C,^ is odd, while all the other fields defined up to this point are even. (That is, all the 
fermions in Table I, the fermion and the Higgs helds are even under P.) The P 

prevents an explicit (or bare) mass term for C,, C,^. But such a mass can arise dynamically as follows. 
Suppose that there are terms C'^C'S', and S, where S' is a singlet scalar and 'h and 

are fermions that transform under a force that confines at a scale A and acquire a condensate 
(xj/ij/) ~ A^. Assume that S and are odd under P while T is even, so as to allow these couplings. 
Then {S) r\j /MQjjrp^ which can be of order 1 GeV if the confinement scale A ~ 10^^ GeV. 

The mechanism just described for making the mass of C, small may seem somewhat contrived, but 
it is “natural” in the technical sense m, i. e. in the sense that no coefficient of a renormalizable term 
in the Lagrangian has to be “fine-tuned” to be very small. Perhaps a somewhat more elegant way 
could be found to make small. But this example does serve to illustrate the general fact that the 
mass of a fermion can be much smaller than the breaking scale of the symmetry that “protects” it 
(here the symmetry P) , without violating the condition of “technical naturalness”. Another example 
of this (which in fact is closely analogous to the mechanism just described) is the type II see-saw 
mechanism for neutrino masses, where neutrino masses are protected by the electroweak SU{2)l, but 
the neutrino masses that result from the breaking of SU{2)l at the scale Fermi scale v are only of 
order /Mgut ^ v. 

With all the particles of the model having the R charges and K and P parities given above, it is 
easy to see that the lowest dimension operator that connects the light dark matter particles C, and 
to the other particles of the low energy theory (i.e. those with mass less than or of order M*) is the 
dim-5 operator which is suppressed by a large {0{Mp^) or 0{Mgut)) mass. That is 

important, as seen in the previous sections, as otherwise there would be a danger that interactions 
of C and with other matter would produce a thermal abundance of them and and lead to a large 
symmetric component of dark matter. 

One can see that the quantum number X serves to stabilize not only the dark matter, but also 
the heavy vectorlike 5-1-5 fermions. This is easiest to see in the case (which we shall henceforth 
assume) that all the vectorlike (VL) fermions, which have odd values of X, are heavier than all the 
DM fermions, which have even values of X. Then the lightest particle with odd X will be absolutely 
stable, and the lightest particle with X = ±2 will also be stable. That X conservation stabilizes both 
the dark matter and the vectorlike fermions seems to be a generic feature of the kind of model we are 
discussing if it is unified in a simple group. If X remained exactly unbroken, there would be residual 
heavy quarks and leptons from these vectorlike multiplets that would come to dominate the energy 
density of the universe. It is thus necessary that some small spontaneous breaking of X in the model. 
It must be large enough for the decay of the heavy vectorlike fermions to be sufficiently fast not to 
cause cosmological problems, but small enough that (a) it does not cause too rapid decay of the heavy 
dark matter Xj (be. dominating over the decays caused by the dim-5 operator), and (b) X-violating 
processes do not wipe out any asymmetries in dark matter. 

Fortunately, there are simple ways that such X violation can be introduced into the model that 
satisfy these conditions. One way is to suppose that there are X-violating VEVs of the Higgs multiplets 

and that are of order e <C 1 times the X-conserving VEVs: (h'j) ~ ev^ ~ and 

(of)) ~ eM*. Then AX = ±1 effects are suppressed by 0(e). This would have the following effects. 

(1) There would be an 0(e) mixing of the SM down quarks d = ^[“^1 with the heavy vectorlike down 
quarks D = so that D could decay into u + W- with an amplitude suppressed by a factor of e. 

(2) There would be similar mixing of the SM charged leptons with the heavy vectorlike charged 
leptons L~, allowing the latter also to decay weakly with an amplitude suppressed by a factor of e. 


(3) There would be mixing among all the neutral fermions, namely the SM neutrinos (i/), the heavy 
vectorlike neutrinos (N), and the heavy dark matter particles (x)- This would allow weak decay of 
the heavy dark matter particles into i~ However, this amplitude would be suppressed by order 

e^. The reason is that the diagonalization of the mass matrices of the charged and neutral fermions 
turns out to give electroweak lepton doublets of the following form: 


/ u + 0{e)N + 0{e^)x \ / N + 0{e)u + 0{e)x 

r + 0 {e)L- j ’ V 


(19) 


The powers of e are directly traceable to the number of units by which X is violated. Because a decay 
of a heavy dark matter particle into Standard Model particles violates X by two units, its amplitude 
is suppressed by e^. 


The condition that the X-violating decay of the heavy dark matter be slow compared to that 
caused by dim-5 operator is that <C M^/MQjjrp => <C M^/Mpi, which is also the condition 

that X-violating scattering be slow enough not to wipe out the dark matter asymmetries. But that 
condition still allows the weak decays of the heavy vectorlike quarks and leptons to happen very 
rapidly; indeed with a lifetime short compared to the Hubble time when T ~ M*. (These issues will 
be analyzed in more detail in a future paper m-) 


It should be noted that despite the spontaneous violation of X, the C, a-i'e extremely stable on 
the time scale of the present age of the universe, since the lowest-dimension operators that couple 
C , to the light particles of the model have dimension 5 and are therefore suppressed by 0{1/Mgut) 
or 0{1/Mp(), and also because their decay amplitude to Standard Model particles would be further 
suppressed by two powers of e. For example, if their decay amplitudes are of order e^/McuT, they 

have a lifetime of ~ ’ which for ttiq ~ 5 GeV and <C M^/Mp£ gives S> 10^® yrs. 


In the SU{7) model we have presented, both the dark matter particles and heavy vectorlike quarks 
and leptons carry non-zero X. (This is a generic feature of such models when unified in a simple group.) 
As one can see from Table I, one can write X = 2Np)M + TVvl, where NpM is dark matter number 
(the number of dark matter particles minus the number of anti-dark matter particles) and Nyi is 
the vectorlike particle number. The contribution of the dark matter fermions and vectorlike fermions 
to XX in S'f7(2)* sphaleron processes cancel, as can be seen from Table I. In other words, SU{2)^ 
sphaleron processes violate Np,M and Nyp, but do not violate X. (This is also a generic feature 
of models of the type we have been discussing when unified in a simple group.) Thus, whatever 
primordial asymmetry is generated after reheating must have X = 2Ni)m + Xlyp ^ 0, as otherwise 
SU{2)^ sphalerons will wipe out the asymmetries in both NpM and Nyp. (This is analogous to the 
fact that electroweak sphalerons do not violate B — L, so that unless a primordial asymmetry in B — L 
exists electroweak sphalerons will wipe out both B and L.) 


4 Conclusions 

We have presented a scenario in which a new non-abelian gauge interaction, whose group we called G*, 
is responsible for both co-generating a dark matter asymmetry and pre-annihilating the symmetric 
component of the dark matter. The co-generation of dark matter with ordinary matter is done by the 
sphaleron processes of G*, while scattering mediated by the exchange of G* gauge bosons accomplishes 
the pre-annihilation of dark matter with anti-dark matter, so as to leave almost purely asymmetric 


dark matter. 


We showed that sufficient pre-annihilation in this scenario requires the scale of G* breaking M* to 
be below about 3 x 10® GeV. The pre-annihilation mechanism requires that there be heavy unstable 
dark matter particles, an asymmetry in which is initially generated and which later decay to a light 
stable form of dark matter, which inherits this asymmetry. (These heavy unstable dark matter parti¬ 
cles, which we call Xi have mass rriy^ ~ M*. The light stable dark matter particles, which we call C, 
should have mass of order 5 GeV to obtain ^dm ~ We have shown that these x decays must 

have lifetimes that fall within a certain relatively narrow range. Shorter lifetimes would cause these 
decays to regenerate the symmetric component of dark matter, while longer lifetimes would result in 
the stable dark matter being too hot. The range for T^ec (the temperature when the y particles decay) 
is roughly between 2 x 10“^ and 10“® times ruy^. We have shown that the x particles can have decay 
lifetimes in the required range, if they come from dim-5 operators suppressed by 1/Mgut or l/Mpg. 

Finally, we have shown that this scenario can be implemented within the context of grand unified 
models based on groups of rank > 6. We presented a fairly simple SU{7) model as an example. This 
model illustrates how several features of the scenario can arise naturally within grand unification; (1) 
Standard-Model-singlet fermions that can play the role of the heavy dark matter are contained in the 
fermion multiplets of the large unified group, thus unifying ordinary quarks and leptons with dark 
matter; (2) global symmetries that stabilize the dark matter can arise as accidental global symmetries, 
analogously to how B — L arises as an accidental symmetry in minimal SU{5); and (3) new (non¬ 
standard Model) non-abelian gauge symmetries are left after the breaking of the grand unified group, 
which can play the role of G*. 

There are several model-building issues that deserve further study. One has to do with the breaking 
of the global symmetry X that stabilizes dark matter. We have found that in fully gauge unified models 
that realize the scenario we discuss here, the global symmetry X not only stabilizes dark matter but 
also stabilizes certain heavy vectorlike fermions that generically exist in such models. X must therefore 
be slightly broken in such a way that the heavy vectorlike particles decay rapidly, while keeping the 
dark matter sufficiently stable. We showed one way to do this, but it would interesting to see if there 
are other better ways. Similarly, the way that the light dark matter particles C were given small mass 
in a “technically natural” way could possibly be improved upon. 


Acknowledgements 


S.M.B. was supported in part by the Department of Energy (DE-FG02-12ER41808). R.J.S. was 
supported in part by the Department of Energy (de-sc0011981). 


References 

[1] S.M. Barr, Phys. Rev. D85, 013001 (2012). 

[2] S.M. Barr and Heng-Yu Ghen, JHEP 1310, 129 (2013). 

[3] V. Kuzmin, V. Rubakov, and M. Shaposhnikov, Phys. Lett. B155, 36 (1985); P. Arnold and L. 
McLerran, Phys. Rev. D36, 581 (1987); ibid. D37, 1020 (1988). 

[4] S.M. Barr, R.S. Chivukula, and E. Farhi, Phys. Lett. B241, 387 (1990). 


[5] S.M. Barr, Phys. Rev. D44, 3062 (1991); D.B. Kaplan, Phys. Rev. Lett. 68, 741 (1992); J. 
Shelton and K.M. Zurek, Phys. Rev. D82, 123512 (2010); N. Haba and S. Matsumoto, Prog. 
Theor. Phys. 125, 1311 (2011); M.T. Frandsen, S. Sarkar, and K. Schmidt-Hoberg, Phys. Rev. 
D84, 051703 (2011). 

[6] M. Blennow, B. Dasgupta, E. Fernandez-Martinez, and N. Rius, JPIEP 1103:014 (2011). 

[7] V.A. Kuzmin, Phys. Part. Nucl. 29, 257 (1998); N. Cosme, L. Lopez-Honorez, and M.H.G. 
Tytgat, Phys. Rev. D72, 043505 (2005); D.E. Kaplan, M.A. Luty, and K.M. Zurek, Phys. Rev. 
D79, 115016 (2009); T. Cohen, D.J. Phalen, A. Pierce, and K.M. Zurek, Phys. Rev. D82, 056001 
(2010); H. An, S.-L. Chen, R.N. Mohapatra, and Y. Zhang, JPIEP 1003:124 (2010); E.J. Chun, 
Phys. Rev. D83, 053004 (2011); J.J. Heckmann and S.J. Rey, JHEP 1106, 120 (2011); N.F. 
Bell, K. Petraki, I.M. Shoemaker, and R.R. Volkas, Phys. Rev. D84, 123505 (2011); C. Cheung 
and K.M. Zurek, Phys. Rev. D84, 035007 (2011); J. March-Russell and M. McCullough, JCAP 
1203, 019 (2012); Y. Cui, L. Randall, and B. Shuve, JPIEP 1204, 075 (2012); C. Arina and N. 
Sahu, Nucl. Phys. B853, 666 (2012); S.M. Barr, Phys. Rev. D84, 013001 (2012). 

[8] S. Nussinov, Phys. Lett. B165, 55 (1985) 

[9] M.L. Craesser, I.M. Shoemaker, and L. Vecchi, JHEP 1110, 110 (2011). 

[10] T. Lin, H.-B. Yu, and K.M. Zurek, Phys. Rev. D85, 063503 (2012). 

[11] K. Petraki and R.R. Volkas, Int. J. Mod. Phys. A28, 1330028 (2013); K.M. Zurek, Phys. Rept. 
537, 91 (2014). 

[12] R. Kitano and I. Low, Phys. Rev. D71, 023510 (2005). 

[13] Zhaofeng Kang and Tianjun Li, JHEP 1210, 150 (2012). 

[14] J. Unwin, JHEP 1306, 090 (2013). 

[15] R.J. Scherrer and M.S. Turner, Phys. Rev. D33, 1585 (1986); D34, 3263 (1986); K. Criest and 
D. Seckel, Nucl. Phys. B283, 681 (1987). 

[16] E. W. Kolb and M. S. Turner, The Early Universe, (New York: Addison-Wesley, 1990). 

[17] S.M. Barr, Heng-Yu Chen, and R.J. Scherrer, in preparation. 

[18] C. ’t Hooft, in Recent Developments in Gauge Theories, ed. C. ’t Hooft et al. (Plenum, New 
York, 1980). 


